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1^ \ Abstract: We study the dependence of the center-projected string tension on both 

'^ ' the lattice size, and the number of gauge copies used for maximal center gauge 

fixing. We show that a recent finding of Bornyakov, Komarov, Polikarpov, and 
r\ \ Veselov ( [hep-lat/0002017D , indicating a substantial breakdown of center dominance 

c^ \ in maximal center gauge, is only obtained for rather small lattice volumes, and is not 

found in numerical simulations on larger lattices. It is shown that center-projected 
Creutz ratios approach the full asymptotic string tension as lattice size increases, and 
that the P- vortex density is consistent, at moderately weak couplings, with 2-loop 
scaling behavior. 
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1. Introduction 

In the last few years there has been renewed interest, within the lattice gauge theory 
community, in the center vortex theory of confinement [1]. The revival of this old 
idea is due to a number of numerical studies, which all indicate that center vortices 
are ubiquitous in the QCD vacuum and give rise to the linear confining potential 
[2-15]. Recently, however, a paper by Bornyakov, Komarov, Polikarpov, and Veselov 
(BKPV) [16] appeared which questions the validity of a procedure, known as center 
projection in maximal center gauge, which was essential in many of these studies. 
The claim is that when large numbers of random gauge copies are used in fixing the 
gauge (in an effort to minimize the Gribov copy problem), the center-projected string 
tension underestimates the full string tension by as much as 30%. Our purpose in the 
present article is to show that, while BKPV have certainly raised an important issue, 
their actual data was strongly affected by the rather small lattice volumes used in 
the numerical simulations. We will show that the conclusions drawn from this data 
are not sustained by simulations on larger lattices, which generally agree with our 
previously reported results [4]. 

In the direct version of maximal center gauge [4] , the procedure is to maximize 

^ = EETr^[f^.(^)] (1-1) 

by an iterative over-relaxation procedure, where Tryi[f/] is the trace of U in the 
adjoint representation. Let R^ denote the value of R after n over-relaxation sweeps. 
When Rn is judged to have converged, e.g. according to a criterion of the form 

Rn^Rn-^^^ (1.2) 

Rn 



then the hnk variables U^{x) are projected onto the nearest center element Z. In 
SU(2) lattice gauge theory, the projection is simply 

Z^{x) = signTr[[/^(x)] (1.3) 

Center-projected Wilson loops, Creutz ratios, etc. are observables computed from 
the center-projected link variables. It was shown in a number of studies (see, in 
particular, ref. [4]), that 

• thin vortex excitations of the projected lattice, known as "P- vortices" , are 
located roughly in the middle of thick center vortices on the unprojected lattice; 

• projected Creutz ratios Xcp{I, I) are close to the asymptotic string tension on 
the unprojected lattice ("center dominance"); 

• the density of P-vortices, from f3 = 2.3 onward, scales according to asymptotic 
freedom; 

• removing center vortices (located via the projected lattice) from unprojected 
lattices also removes confinement and chiral symmetry breaking, and brings 
the topological charge on the lattice to zero [6]. 

Because these numerical results have such strong implications for the QCD con- 
finement mechanism, it is important that the center projection procedure be exam- 
ined critically. The first obvious question — why should this procedure work at all? 
— was addressed in ref. [17]. There it was shown that in the absence of Gribov 
copies (i.e. if the gauge can be fixed to a global maximum of R), then center projec- 
tion in maximal center gauge will always locate a thin vortex inserted anywhere on 
the lattice. This was dubbed the "vortex-finding property" of maximal center gauge, 
and is certainly a necessary condition for its success. However, the vortices found 
in the QCD vacuum are not thin vortices, but are necessarily of finite thickness in 
physical units. Moreover, maximal center gauge is plagued with Gribov copies, since 
the over-relaxation scheme converges only to a local maximum of R, which will be 
slightly different for every gauge copy of a given lattice configuration. 

In view of the Gribov copy problem, and the finite thickness of vortices, one must 
rely on empirical checks of vortex-finding via center projection. In this article we will 
study the sensitivity of center-projected Creutz ratios and the P-vortex density with 
respect to: (i) the number N^opy of random gauge copies used for maximal center 
gauge fixing; (ii) the lattice size; and (iii) the convergence parameter 6 in eq. (|1.2| ). 
We will here only discuss center projection in direct maximal center gauge, since this 
is the case treated in ref. [16], and it is also the center gauge with which we have the 
most experience. It should be noted, however, that alternatives to maximal center 
gauge do exist; these include Laplacian center gauge [7] (which is free of the Gribov 
copy problem), as well as two other recent proposals [18,19]. 



2. Ncopy and Lattice Size Dependence 

We begin with the Ncopy dependence, whose importance was recently emphasized 
by BKPV [16]. As noted above, when the over-relaxation gauge-fixing procedure 
is apphed to different gauge copies of a given lattice configuration, slightly different 
values of R are obtained. One way to minimize the gauge copy dependence is to carry 
out the over-relaxation procedure on a number Ncopy of random gauge copies, perform 
center projection on the copy with the largest value of R, and evaluate observables. 
Data obtained in this way, over a range of Ncopy values, can then be extrapolated to 
the Ncopy -^ oo limit. ^ In the original simulations of ref. [4], only three gauge copies 
were used, and there was no attempt to extrapolate to Ncopy -^ oo. Projected Creutz 
ratios Xcp{I, I) were found to be close to the asymptotic string tensions reported by 
Bali et al. [21], for all / > 2. BKPV, however, calculate Creutz ratios in the range 
Ncopy £ [1, 20], and extrapolate to Ncopy ^ oo by fitting their data to the functional 
form 

xZ-^^{I,I)=Xcp{IJ) + ^^ (2.1) 

y^^copy 

The result reported by BKPV is that projected string tensions, at (3 = 2.4, 2.5, 
underestimate the full string tension by about 20% at Ncopy = 20, and by as much 
as 30% in the extrapolation to Ncopy — ^ oo. This result suggests that the center 
dominance previously reported in maximal center gauge was a numerical accident, a 
result of using too few gauge copies. 

However, apart from the number of gauge copies, there is one other notable 
difference between the BKPV simulations and previous work. This is the matter 
of lattice size. Data reported in ref. [4] was obtained at /3 = 2.3 and (3 = 2.4 
on 16^ lattices, and at (3 = 2.5 on a 22*^ lattice. BKPV, on the other hand, used 
only 12^ lattices at /5 = 2.3 and (3 = 2.4, and lattice size 16^ at (3 = 2.5. This 
raises the question of whether the BKPV results, obtained on the smaller lattices, 
were seriously contaminated by finite-size effects. To find out, we have repeated the 
center-projection calculation at (3 = 2.3 and f3 = 2.5 on a variety of lattice sizes, for 
Ncopy £ [1, 20]. For the convergence parameter in eq. (|1.2|) , we have used 6 = 2x 10"'''. 

In Fig. I we display results for the Creutz ratio x^''°™(4, 4) vs. Ncopy at /3 = 2.5, 
for lattice sizes ranging from 8^ to 28^. Two features of this data are immediately 
apparent. First, there is indeed a slow downward trend in the Creutz ratio as Ncopy 



^P-vortex positions extracted from different Gribov copies are closely correlated, as one would 
expect if they locate physical objects [4]. It is possible, however, by choosing a very special starting 
configuration with links gauge-fixed to have mainly positive trace, to destroy the vortex-finding 
property of center projection even for thin vortices [17, 20]. This is one aspect of the Gribov copy 
problem in maximal center gauge. Since the vortex-finding property is essential to center projection, 
it is evident that any unusual variation of the iterative gauge-fixing procedure, which destroys this 
property, must be avoided. 
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Figure 1: Center projected Creutz ratios Xcp°''^(4, 4) vs. Ncopy at various lattice volumes, 
for /3 = 2.5. The solid line in this figure, and in Figs. |2[S indicates the asymptotic string 
tension extracted by standard methods on unprojected lattices, reported by Bali et al. [21]. 
Dashed lines indicate the errorbar in this asymptotic string tension. 



increases, as noted by BKPV, but this effect is much more pronounced on smaller 
lattices than on larger lattices. Second, although Creutz ratios on the smaller lattices 
grossly underestimate the full string tension, the data appears to steadily increase 
towards the full asymptotic string tension, reported by Bah et al. [21], as the lattice 
size increases. These trends in the data are by no means unique to the particular 
Creutz ratio Xcp(4, 4) at /5 = 2.5, but are typical of all of our results. For completeness 
we display, in Figs. ^ and ^ some other projected Creutz ratios X^''°™(-^) ^) fo^' ^ ^^ 
the range / = 2 — 5, at couplings P = 2.3 and /5 = 2.5. In Figs. ^^ solid (dashed) 
lines indicate the value (errorbar) of the asymptotic string tension on the unprojected 
lattice, reported in ref. [21]. 

In Fig. ||we show the projected Creutz ratios X^''°™(/, /) for N^opy = 5, 10, 15, 20 
on the largest lattices we have used: 20"^ at /9 = 2.3 and /? = 2.4, and 28^ at /5 = 2.5. 
We also show the values of these Creutz ratios extrapolated to the Ncopy — ^ oo limit, 
using the fitting function ( p.l| ) suggested by BKPV. As usual in maximal center 
gauge, all the Xcp{I, I) for / > 2 are close to the asymptotic string tension, and these 
latest results are not far from our earlier results reported in ref. [4] . As another way 
of showing lattice size dependence, we take the average of the projected (A^, 
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Figure 2: Center projected Creutz ratios Xcp°^^{I,I) vs. Ncopy and lattice volume, at 
/? = 2.3 and / = 2 - 5. 



Creutz ratios Xcp{I, I) in the range 1 = 2 — 5 

1 5 

Xav = T Z^ Xcpl-' 5 -'J 
^ 7=2 



(2.2) 



Figure | shows the lattice size dependence of Xav Note the approach of Xav to the 
full asymptotic string tension {crsu{2) = 0.035 at /5 = 2.5) as lattice size increases. 

2.1 Gauge-Fixing Convergence Criterion 

In addition to lattice volume and Ncopy dependence, it is also worthwhile to check 
that the numerical results are stable when the gauge-fixing convergence criterion is 
strengthened, i.e. when the constant S in eq. (|1.2|) is reduced. When 6 is chosen too 
large, the center projected Creutz ratios come out significantly too high. This is il- 
lustrated in Fig. 1^, which shows results for projected Creutz ratios with convergence 
criteria 6 = lO^^ 10~^ 10-^ 2 x 10"^ at /3 = 2.5 (24^ lattice). The weakest conver- 
gence criterion, corresponding to 6 = 10~^, is clearly insufficient for accurate results, 
but Creutz ratios obtained with the two smallest values of S are fairly consistent, 
indicating that these numbers are not far from the 5 — > limit. 
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Figure 3: Center projected Creutz ratios Xcp°^^{I,I) vs. N^opy and lattice volume, at 
/? = 2.5 and / = 2 - 5. 



3. Vortex Density 

Next we consider the scaling of the vortex density. The lattice P-vortex density p 
is the total number of P-vortex plaquettes (i.e. plaquettes on the projected lattice 
with ZZZZ = —1), divided by the total number of plaquettes on the lattice. This 
quantity is proportional to the average area taken up by P-vortices per unit lattice 
volume, and is determined from the center-projected plaquette expectation value 



p=^(l-W,,[l,l]) 



(3.1) 



If this quantity scales as predicted by asymptotic freedom, then we would have 
P = Paf, where 



Paf 



P 
6A2 



Fi(3) 



(3.2) 



Here p is the vortex density (average vortex area per unit volume) in physical units 
of inverse area, and 
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Figure 4: Center-projected Creutz ratios Xcp{I, I) at A^copy 
lated to Ncopy — > oo, on the largest lattices (20'* lattices at /3 = 
used in our simulations. 



= 5, 10, 15, 20, and extrapo- 
2.3,2.4, and 28^ at /3 = 2.5) 



In Fig. we plot the lattice vortex density rescaled by the asymptotic freedom 
expression 



P = 



P 



Fm 



(3.4) 



which should be constant in the large (3 limit, if p scales according to asymptotic 
freedom. There appears to be good evidence for this kind of scaling, already for 
(3 > 2.2, in agreement with previous results [4, 10]. In Fig. ^ and Table |I|, the vortex 
densities at /5 = 2.3, 2.4, 2.5 are taken from the largest lattices (20^, 20^ and 28^, 
respectively), and extrapolated to the Ncopy -^ oo limit by a fit to 
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The vortex densities at other values of /5 (with Ncopy = 3) are just taken from our 
previous work. For comparison, in Table 0, we display our values for p, and the 
values of the rescaled asymptotic string tension 
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Figure 5: Average of projected, Ncopy -^ c« extrapolated Creutz ratios Xcpil,!), in tiie 
range 7 = 2 — 5. 
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Figure 6: Effect of varying the convergence criterion 6 on projected Creutz ratios (ex- 
trapolated to Ncopy -^ oo) at /? = 2.5 on a 24^ lattice. 



where crsui2) is the string tension (in lattice units) on the unprojected lattice, at the 
given f3 value. It is interesting to note that the scaling of vortex density p, in the 
range f3 = 2.3 — 2.5, is substantially better than the the scaling of the full asymptotic 
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Table 1: Rescaled vortex density p and SU(2) string tension a. These quantities should 
be constant in the scaling limit. 

string tension asu{2) in this range. 
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Figure 7: Rescaled vortex density p = p/F{P), where p is the measured vortex density 
and F{P) is the asymptotic freedom expression in eq. (|3.3|) . Vortex density scales Up is 
constant. 



4. Vortex Thickness 



Why are such relatively large lattices required, in order to avoid large finite size effects 
in center projected Creutz ratios and vortex densities? We believe that the relevant 
length scale here is associated with the vortex thickness. Center vortices are surface- 
like objects, in D = 4 dimensions, which have a finite thickness in physical units. 
The thick vortex surface (or "core") bounds a Dirac 3- volume, which represents the 
region of discontinuity of a singular gauge transformation associated with the vortex. 
In ref. [17] we have explained the vortex- finding property of maximal center gauge 
in terms of the global properties of this singular gauge transformation. On a small 



lattice, however, with an extension comparable to the vortex thickness, these global 
aspects of the vortex field may be almost absent, and the minimal Dirac 3- volume 
could be quite small. In that case, the argument of ref. [17] breaks down. For this 
reason, we expect center projection to be less effective at finding vortices on small 
lattices, leading to underestimates of both the vortex density and the projected 
Creutz ratios. 

Assuming there is some truth in this explanation, center projection can only be 
accurate for lattices whose extension is large compared to the vortex thickness. There 
are three independent ways of estimating the thickness of center vortices, which can 
be deduced from either 

1. the ratio of "vortex-limited" Wilson loops [4]; 

2. the vortex free energy as a function of lattice size [14]; 

3. the adjoint string-breaking length [22]. 

Vortex-limited Wilson loops are defined in the following way: Wn{C) is a Wilson loop 
evaluated on a sub-ensemble of unprojected configurations, selected so that precisely 
n P-vortices, in the corresponding center-projected configurations, pierce the minimal 
area of the loop. We can further make the restriction, for Wi{C), that the negative 
P- vortex plaquette hes at (or touches) the center point of the rectangular loop. It is 
then expected that 

^5^ ^ -1 (4.1) 

in the limit where the vortex core is entirely contained within the loop (cf. ref. [4] 
for a more extended discussion). In Fig. H we show the data for Wi/Wq vs. loop area 
at /5 = 2.3, taken from our previous work in ref. [4]. Judging from this figure, the 
vortex appears to almost fit inside a 5 x 5 loop, which leads to a rough estimate of 
the vortex radius, as it pierces a plane, of about 3 lattice spacings. At /5 = 2.3 we 
have aa^ = 0.135, and taking a = 5 fm~ , the lattice spacing is a = 0.164 fm. A 
diameter of 6 lattice spacings then corresponds to a vortex thickness of about one 
fermi. 

A second estimate is obtained from the recent calculation of vortex free energy 
vs. lattice size, carried out numerically by Kovacs and Tomboulis [14]. The vortex 
free energy is close to zero when the lattice extension is greater than the vortex 
thickness, and this again gives an estimate for the vortex thickness of a little over 
one fermi. Finally, if confinement is due to center vortices, then an i? x T Wilson loop 
in the adjoint representation must change from a (Casimir scaling) area-law falloff 
to a (color-screening) perimeter-law falloff for charge separation R greater than the 
vortex thickness [23]. The adjoint string-breaking distance has been measured, by de 
Forcrand and Philipsen [6], to be 1.25 fm, and this distance provides us with a third 
estimate of the vortex thickness, which is roughly consistent with the other two. 
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Figure 8: Ratio of the 1-Vortex to the 0- Vortex Wilson loops Wi{C)/Wo{0), vs. loop area 
at /? = 2.3. 

At /3 = 2.5, one fermi corresponds to 12 lattice spacings. The lattice used by 
BKPV at /? = 2.5 was only 16 lattice spacings across, and this may simply be 
inadequate for center projection to reliably identify vortices, in view of the above 
estimates for the vortex thickness. In fact, the P- vortex density is quite low on small 
lattices, increasing sharply up to L = 16 lattice spacings at /5 = 2.5, where it begins 
to level off. This is illustrated in Fig. ^; the data points are the extrapolated values 
for p at A*", 
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In order to estimate the average distance between vortices, we need to know the 
vortex density. The P-vortex density, discussed in the previous section, is in fact an 
overestimate of the actual center vortex density [10,12]. The reason is that P- vortices 
fluctuate within the thick vortex core, and P-vortex locations, while correlated among 
random Gribov copies, do vary somewhat from one random copy to another [4]. Thus, 
although P- vortices are certainly near the middle of center vortices, they are unlikely 
to be exactly in the middle. A more accurate estimate of the center vortex density is 
arrived at by either "smoothing" the P-vortex surfaces, as explained in ref. [12], or 
else by taking the naive estimate of the number of vortices piercing a plane (per unit 
area in lattice units), extracted directly from the string tension (see, e.g., ref. [23]) 



/=2(1 



(4.2) 



The two estimates agree fairly well [12], and for jS = 2.3 we find / = 0.063. This 
implies an average distance of /^^^^ ~ 4 lattice spacings between the centers of 
vortex cores piercing a plane. Since we have already estimated the vortex thickness 
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Figure 9: Vortex density p vs. tiie extension of L of the iiypercubic lattice (volume L^) 
at 13 = 2.5. 

aX f3 = 2.3 to be about 6 lattice spacings, its clear that there must be a substantial 
overlap between vortex cores, even on a very large lattice. There is nothing in 
principle wrong with that; vortex cores are not impenetrable objects, and their long- 
range effects are associated with Dirac 3- volumes, rather than the detailed structure 
of the core. Our findings for vortex thickness and separation simply indicate, in 
accordance with some old ideas of the Copenhagen group (Nielsen, Olesen, and 
Ambj0rn in ref. [1]), that the QCD vacuum is more like a liquid of vortices than 
a dilute gas. Perhaps this is natural for objects which are both condensed and 
disordered. 

5. Conclusions 

We have found that center-projected lattices are more sensitive to finite size effects 
than are unprojected lattices, and precision results for center-projected Creutz ratios 
require lattice sizes which are large compared to one fermi. If accurate numbers 
(rather than just qualitative results) are required, then center-projected data from 
numerical simulations must be tested for convergence with respect to: (i) increasing 
the lattice size; (ii) increasing the number Ncopy of gauge copies used for center gauge 
fixing, and (iii) strengthening the gauge-fixing convergence criterion. 

Our work was stimulated by the recent findings of BKPV [16], who have demon- 
strated a significant A^copy-dependence for center-projected Creutz ratios obtained on 
relatively small lattice volumes. In the large Ncopy hmit, a large deviation {^ 30%) 
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was found between the string tension on unprojected and projected lattices. Our 
result in the present article is that this Ncopy dependence is greatly reduced as lattice 
size increases, and center projected Creutz ratios Xcp{I, I) appear to converge to val- 
ues which are quite close to the asymptotic string tension obtained on unprojected 
lattices. 

It is difficult to say whether the full and center-projected string tensions are in 
precise agreement in the infinite volume limit. To check such agreement, what is 
really required is an extrapolation to a triple limit: volume — *> oo, Ncopy -^ oo, and 
S —^ 0; our data is not yet adequate to extrapolate to this triple hmit systematically. 
But we emphasize again that on the largest lattices we have used, and with results 
extrapolated (following BKPV) to the Ncopy — ^ oo limit, the center projected Creutz 
ratios lie quite close to the asymptotic string tension, as illustrated in Fig. | above. 
Moreover, on the same large lattices and Ncopy —^ oo extrapolation, there is good 
evidence for asymptotic scaling of the vortex density at couplings P > 2.2. The 
scaling of the vortex density, at the /? values studied, is in fact substantially better 
than the scaling of string tension itself. 

Center projection in maximal center gauge is intended as a method for locating 
center vortices, on unprojected lattices, from the position of P-vortices on projected 
lattices. The correspondence of P-vortices with physical objects depends on 

1. Scaling of the P- vortex density; 

2. Correlation of P-vortex locations on the projected lattice with gauge-invariant 
observables on the unprojected lattice. 

In particular, it is the correlation of P-vortices with unprojected Wilson loops, and 
specifically the ratios of "vortex-limited" Wilson loops Wn{C)/Wo{C) on the unpro- 
jected lattice, which indicate that P-vortices correspond to center vortices (rather 
than to some other type of object) on the unprojected lattice [5]. Moreover, if P- 
vortices locate center vortices, and if center vortices produce the full asymptotic 
string tension, then we must find 

3. Center dominance; i.e. the equality of the string tension on projected and 
unprojected lattices. 

The numerical simulations reported here were concerned with the scaling of the P- 
vortex density, and with center dominance. From the results of those simulations, we 
conclude that center projection in maximal center gauge remains a viable method of 
locating center vortices on full, unprojected lattices, and that these vortices probably 
account for the entire asymptotic string tension. 



13 



Acknowledgments 

Our research is supported in part by Fonds zur Forderung der Wissenschaftlichen 
Forschung P13997-PHY (M.F.), the U.S. Department of Energy under Grant No. 
DE-FG03-92ER40711 (J.G.), the "Action Austria-Slovakia: Cooperation in Science 
and Education" (Project No. 30sl2) and the Slovak Grant Agency for Science, Grant 
No. 2/7119/2000 (S.O.). J.G. also acknowledges the support of MaPhySto, Centre 
for Mathematical Physics and Stochastics, funded by the Danish National Research 
Foundation. 

References 

[1] G. 't Hooft, Nucl. Phys. B138 (1978) 1; 

G. Mack, in Recent Developments in Gauge Theories, edited by G. 't Hooft et al. 

(Plenum, New York, 1980); 

H. B. Nielsen and P. Olesen, Nucl. Phys. B160 (1979) 380; 

J. Ambj0rn and P. Olesen, Nucl. Phys. B170 (1980) 60; 265; 

P. Vinciarelli, Phys. Lett. B78 (1978); 

T. Yoneya, Nucl. Phys. B144 (1978) 195; 

J. M. Cornwall, Nucl. Phys. B157 (1979) 392; 

R. P. Feynman, Nucl. Phys. B188 (1981) 479. 



[9 

[lo: 



L. Del Debbio, M. Faber, J. Greensite, and S. Olejnik, Phys. Rev. D55 (1997) 2298, 
|hep-lat/9610005| . 

L. Del Debbio, M. Faber, J. Greensite, and S. Olejnik, in New Developments in Quan- 
tum Field Theory, ed. Poul Henrik Damgaard and Jerzy Jurkiewicz (Plenum Press, 
New York-London, 1998) 47, piep-lat/9708023. 



L. Del Debbio, M. Faber, J. Giedt, J. Greensite, and S. Olejnik, Phys. Rev. D58 (1998) 



094501, |hep-lat/9801027 . 



M. Faber, J. Greensite, and S. Olejnik, JHEP 01 (1999) 008, [hep-lat/9810008 



Ph. de Forcrand and M. D'Elia, Phys. Rev. Lett. 82 (1999) 4582, |hep-lat/990l"020| . 



C. Alexandrou, M. D'Elia, and Ph. de Forcrand, Nucl. Phys. Proc. Suppl. 83-84 (2000) 
437, |hep-lat/990702^ ; and Nucl. Phys. A663-664 (2000) 1031, |hep-lat/9909005| . 



K. Langfeld, H. Reinhardt, and O. Tennert, Phys. Lett. B419 (1998) 317, ^epj 
lat/9710068|. 



M. Engelhardt, K. Langfeld, H. Reinhardt, and O. Tennert, Phys. Rev. D61 (2000) 



054504, |hep-lat/9904004 



M. Engelhardt, K. Langfeld, H. Reinhardt, and O. Tennert, Phys. Lett. B452 (1999) 
301, |hep-lat/9805002. 



14 



[11] 

[12] 



J. Gattnar, K. Langfeld, A. Schafke, and H. Reinhardt, |hep-lat/0005016| . 



R. Bertie, M. Faber, J. Greensite, and S. Olejnik, JHEP 03 (1999) 019, ^ 
lat/9903023| . 



[13] 

[14] 

[15] 
[16] 

[17] 



B. Bakker, A. Veselov, and M. Zubkov, Phys. Lett. B471 (1999) 214, [hep-lat/99020To| ; 
M. Chernodub, M. Polikarpov, A. Veselov, and M. Zubkov, Nucl. Phys. Proc. Suppl. 
73 (1999) 575, |hep-lat/9809T58. 



T. Kovacs and E. Toniboulis, hep-lat /0002004 ; 

E. Toniboulis, talk at ConHnement 2000, Osaka, Japan. 



A. Hart, B. Lucini, Z. Schram, and M. Teper, JHEP 06 (2000) 040, |liep-lat/000501C| . 
V. Bornyakov, D. Komarov, M. Polikarpov, and A. Veselov, JETP Lett. 71 (2000) 



231, |hep-lat/0002017 



M. Faber, J. Greensite, S. Olejnik, and D. Yamada, JHEP 12 (1999) 012, ^ 



lat/9910033 



[18] 
[19] 

[20] 
[21] 
[22] 
[23] 



B. Bakker, A. Veselov, and M. Zubkov, Hiep-lat/000702^ . 

K. Langfeld, talk at Quark ConBnenient and the Hadron Spectrum IV, Vienna, July 
3-8, 2000. 



T. Kovacs and E. Toniboulis, Phys. Lett. B463 (1999) 104, |hep-lat/990502£ . 



G. Bah, C. Schlichter, and K. Schilling, Phys. Rev. D51 (1995) 5165, [hep-lat /9409005 . 

P. de Forcrand and O. Phihpsen, Phys. Lett. B475 (2000) 280, |hep-lat/9912050 . 

M. Faber, J. Greensite, and S. Olejnik, Phys. Rev. D57 (1998) 2603, |hep-lat/9710039| ; 



M. Faber, J. Greensite, and S. Olejnik, Acta. Phys. Slov. 49 (1999) 177, ^epj 
lat/9807008i 



15 



